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Asymptotic  solutions  for  the  nonlinear,  nonhomogeneous,  Korteweg- 
deVries  (KdV)  partial  differential  equation  <pde)  with  slowly  varying 
coefficients  are  not  in  general  uniformly  valid.  A  uniform  asymptotic 
expansion  is  obtained  by  finding  separate  expansions  for  different 
regions  and  matching.  A  KdV  solitary  wave  propagating  in  slowly  varying 
media  is  examined.  Quasi-stationarity  for  the  core  reduces  the  problem 
to  solving  ordinary  differential  equations  for  that  region.  However,  in 
the  leading  tail  region,  hyperbolic  pde's  must  be  solved  to  determine 
the  amplitude  and  phase.  The  method  of  characteristics  predicts  triple 
valuedness  after  a  caustic  (penumbral  or  cusped)  develops.  Singular 
perturbation  methods  show  the  solution  near  first  focusing  satisfies  the 
diffusion  equation  and  involves  either  an  incomplete  Airy-type  integral 
or  an  exponential  integral  similar  to  the  Pearcey  integral.  Laplace's 
method  shows  that  the  critical  points  of  the  exponential  phase  satisfy 
the  fundamental  folding  equation.  A  linear  multi-phase  solution  is 
determined  which  does  not  become  triple  valued  (break).  Instead,  a  wave 
number  shock  develops,  which  separates  two  different  solitary  wave 
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CHAPTER  I 


INTRODUCTION 


1 .1  History 

The  objective  of  this  thesis  is  to  add  to  our  understanding  of  a 
physical  phenomenon  which  was  first  scientifically  observed  over  150 
years  ago  by  0.  Scott  Russell.  The  incident  in  1834  is  best  described 
by  Russell  himself: 

"I  was  observing  the  motion  of  a  boat  which  was  rapidly  drawn 
along  a  narrow  channel  by  a  pair  of  horses,  when  the  boat 
suddenly  stopped--not  so  the  mass  of  water  in  the  channel 
which  it  had  put  in  motion;  it  accumulated  round  the  prow 
of  the  vessel  in  a  state  of  violent  agitation,  then 
suddenly  leaving  it  behind,  rolled  forward  with  great 
velocity,  assuming  the  form  of  a  large  solitary  elevation, 
a  rounded,  smooth  and  well-defined  head  of  water,  which 
continued  its  course  along  the  channel  apparently  without 
change  of  form  or  diminution  of  speed.  I  followed  it  on 
horseback,  and  overtook  it  still  rolling  on  at  a  rate  of 
some  eight  or  nine  miles  an  hour,  preserving  its  original 
figure  some  thirty  feet  long  and  a  foot  to  a  foot  and  a 
half  in  height."  [30] 

Russell's  life  was  dramatically  changed  by  this  single  event,  and  so  was 
our  understanding  of  wave  motion  to  change  as  the  investigation  began. 

He  found  no  mathematical  theory  available  which  predicted  such  a 
phenomenon  and  proceeded  to  study  it  experimentally.  At  the  time  he 
could  hardly  have  imagined  the  controversy  that  would  arise  between  him 
and  Airy  over  the  existence  of  such  a  permanent  waveform,  much  less  have 
appreciated  the  bountiful  applications  in  science  that  would  arise  and 
the  number  of  great  minds  that  would  apply  themselves  to  broadening  our 
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understanding  of  this  phenomenon!  It  was  nearly  fifty  years  before  the 
controversy  was  put  to  rest  by  such  men  as  Boussinesq,  Rayleigh,  and 
Korteweg  and  deVries.  They  demonstrated  the  existence  of  a  permanent 
solitary  wave  for  nonlinear  partial  differential  equations  of  shallow 
water  theory.  There  was  little  scientific  activity  for  another  half 
century  until  other  physical  applications  were  shown  to  be  governed  by 
similar  mathematics.  Then  a  variety  of  situations  were  shown  to  give 
rise  to  nonlinear  dispersive  partial  differential  equations  (collision- 
free  hydromagnetic  waves  by  Gardner  and  Morikawa  and  a  nonlinear  meson 
field  theory  by  Perring  and  Skyrme).  Fermi,  Pasta  and  Ulam's  study  of 
heat  conductivity  in  solids  then  motivated  Zabusky  and  Kruskal's 
important  numerical  study  of  solitary  waves,  coined  "solitons"  by  them 
because  they  observed  particle-like  interactions.  In  order  to  under¬ 
stand  these  interactions,  a  theoretical  investigation  was  undertaken 
resulting  in  one  of  the  most  significant  developments  in  applied  mathe¬ 
matics  and  mathematical  physics  in  recent  years.  Gardner,  Greene, 

Kruskal  and  Miura  related  inverse  scattering  for  the  Schrodinger  eigen¬ 
value  problem  to  the  nonlinear  Korteweg-deVries  (KdV)  equation.  In  this 
way  they  obtained  solitons  and  multiply-interacting  solitons.  Shortly 
thereafter  Lax  showed  that  other  nonlinear  partial  differential 
equations  could  be  analyzed  in  this  way.  A  series  of  interesting 
physical  problems  was  studied  ( self- induced  transparency  by  McCall  and 
Hahn  collaborated  by  Lamb,  and  the  nonlinear  Schrodinger  equation  by 
Zakharov  and  Shabat,  both  important  examples  since  they  introduced  a 
different  eigenvalue  problem).  Finally,  Ablowitz,  Kaup,  Newell  and 
Segur  showed  how  the  inverse  scattering  transform  (1ST)  could  be  used  to 


solve  a  wide  class  of  nonlinear  evolution  equations.  Thus  Russell's 
discovery  of  the  solitary  wave  has  led  to  much  research  and  no  doubt 
will  lead  to  much  much  more. 

1.2  Motivation 

A  search  of  the  literature  shows  much  activity  and  excitement  in 
the  area  of  solitary  waves  and  solitons  once  Gardner,  Greene,  Kruskal 
and  Miura  [4]  demonstrated  their  method  for  obtaining  exact  solutions  to 
the  nonlinear  KdV  equation.  Their  method  of  direct  and  inverse 
scattering,  which  views  the  eigenfunctions  as  the  transmitted  portion  of 
a  wave  coming  in  from  infinity  and  being  scattered  by  the  initial 
conditions,  was  formalized  by  Ablowitz,  Kaup,  Newell  and  Segur  in  1974 
[2].  This  formalization  demonstrated  permanent  solitary  wave  solutions 
for  a  large  class  of  problems.  The  technique,  however,  requires  that 
the  problem  be  integrable,  which  is  not  the  case  for  many  problems  which 
involve  slow  variations.  It  is  precisely  these  type  problems  which 
arose  in  many  physical  applications.  Researchers  thus  began  using 

I 

numerical  techniques  to  gain  understanding,  and  singular  perturbation 
methods  to  analytically  demonstrate  asymptotic  behavior.  In  the  late 
1970's  several  independent  studies  showed  the  formation  of  a  "shelf" 
behind  a  slowly  varying  KdV  solitary  wave.  Kaup  and  Newell  [14]  and 
Knickerbocker  and  Newell  [19]  used  perturbations  on  the  inverse 
scattering  solution,  while  Ko  and  Kuehl  [16,17],  Karpman  and  Maslov 
[13],  and  Grimshaw  [5]  used  direct  perturbation  methods.  For  the  case 
of  a  KdV  solitcn  propagating  to  the  right  through  a  region  of  slowly 
changing  depth,  Kaup  and  Newell  write: 


"if  a  soliton  propagates  through  a  region  of  decreasing 
depth  to  a  new  constant  level  h-j  at  x  =  x-]  ,  the 
solution  in  the  region  to  the  right  of  x  :  x-|  will 
consist  of  (1)  the  original  soliton,  whose  height  has 
adjusted  in  the  transition  region  so  as  to  be  just  the 
correct  height  for  the  new  depth  plus  (2)  a  number  of 
secondary  solitons  produced  by  the  arrival  at  x  =  x-j  of 
the  elevated  shelf  (which  may  be  treated  approximately  as 
a  rectangular  well  potential  of  width  L  :  (x  -  xg)  and 
average  height  H  *  -L/3n;  the  number  of  solitons  will  be 
proportional  to  /HL)  and  (3)  some  left-over  radiation." 

[14] 

It  is  this  case  of  a  solitary  wave  propagating  on  the  surface  of  water 
of  variable  depth  which  Grimshaw  [5]  considers  using  a  multiple  scale 
approach.  He  shows  that  the  amplitude  and  phase  of  slowly  varying 
solitary  waves  are  determined  by  hyperbolic  partial  differential 
equations.  In  his  paper  he  points  out  for  the  region  in  front  of  the 
wave: 

"A  caustic  develops. . .Beyond  the  caustic  fg  becomes 
multivalued,  and  it  is  tempting  to  conjecture  that  the 
breakdown  is  associated  with  the  formation  of  further 
solitary  waves."  [5] 

It  is  here  that  we  pick  up  the  study.  The  development  of  a  shelf  behind 
a  slowly  varying  solitary  wave  has  been  extensively  studied  (see  above), 
but  the  situation  in  which  characteristics  cross  in  front  of  the  wave 
has  not  been  analytically  explained.  The  primary  impetus  for  the  study 
is  thus  to  explain  analytically  what  the  development  of  a  caustic  in 
front  of  the  wave  means  in  terms  of  the  physical  solution.  Is 
Grimshaw' s  conjecture  correct? 


1.3  Results 

First  we  show  that  the  problem  separates  asymptotically  into  two 
parts.  We  call  the  two  parts  of  the  slowly  varying  solitary  wave  the 


"core”,  or  the  region  where  the  mass  is  concentrated,  and  the  "leading 
tail",  or  exponentially  decaying  region  in  front  of  the  core  with 
respect  to  its  motion.  Our  first  results  came  by  making  some 
simplifying  assumptions.  Looking  first  to  the  leading  tail  we  linearize 
the  problem  and  then  consider  the  special  case  of  constant  coefficients, 
with  slow  variations  introduced  through  initial  conditions  only.  We 
use  the  methods  of  characteristics,  multiple  scales,  matched  asymptotic 
expansions  and  asymptotics  of  exponential  integrals.  The  results  are 
quite  interesting.  Further  solitary  waves  are  not  created  in  the 
leading  tail  as  conjectured  by  Grimshaw.  Instead,  we  find  that  the 
triple  valuedness  predicted  by  the  method  of  characteristics  actually 
leads  to  a  "wave  number  shock",  or  a  relatively  rapid  transition  from 
one  slowly  varying  wave  tail  to  another.  Details  of  these  results  are 
given  in  Haberman  and  Allgaier  [8].  The  idea  of  wave  number  shocks  was 
first  discussed  by  Howard  and  Kopell  [10]  in  their  work  on  slowly 
varying  reaction-diffusion  equations.  Also  Haberman  and  Sun  [9]  used 
similar  asymptotic  calculations  in  their  work  on  slowly  varying 
oscillatory  dispersive  waves. 

The  remaining  chapters  of  this  thesis  show  how  the  properties  of 
the  leading  tail  relate  to  the  structure  of  the  solitary  wave  core. 
Futhermore,  the  most  general  cases  of  arbitrary  slowly  varying  co¬ 
efficients  and  arbitrary  slowly  varying  initial  conditions  for  both  the 
core  and  the  leading  tail  are  solved. 

We  show,  in  fact,  Grimshaw's  conjecture  (see  Section  1.2)  is 
partially  correct,  that  under  the  right  conditions  characteristics  for 
the  leading  tail  will  cross  and  lead  to  the  formation  of  a  new  type  of 
solitary  wave  tail.  We  show  that  the  set  of  different  conditions,  which 
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lead  to  a  wave  number  shock  in  the  leading  tail  of  a  solitary  wave, 
contains  the  particular  conditions  which  Smyth  [2]  demonstrates  lead  to 
the  formation  of  further  solitary  waves  in  the  region  behind  the  core. 

Finally,  we  find  some  very  interesting  facts  concerning  the 
formation  of  caustics.  The  specific  slowly  varying  coefficients  and 
slowly  varying  initial  conditions  are  what  determine  the  type  of  caustic 
which  arises.  In  general,  a  cusped  caustic  (so  named  because  of  its 
shape,  see  Figure  4.1)  occurs  when  the  crossing  of  characteristics  is 
approximated  by  a  complete  cubic  equation.  This  can  occur  both  for 
characteristics  generated  by  the  moving  core  and  for  characteristics 
generated  by  initial  conditions  in  the  leading  tail.  On  the  other  hand, 
a  penumbral  caustic  (here  penumbral  means  partial  illumination  resulting 
from  the  cut-off  of  some  rays  or  characteristics  as  in  an  eclipse) 
occurs  when  the  crossing  of  characteristics  is  approximated  by  an 
incomplete  (or  cut-off)  quadratic  or  cubic  equation.  This  only  occurs 
when  the  characteristic  eminating  from  the  initial  moving  core  catches 
up  to  the  characteristic  from  the  nearest  initial  leading  tail  (see 
Figures  3.1  and  3.3).  For  much  more  detail  on  the  classification  of 
caustics  we  refer  the  reader  to  the  work  by  Kravtsov  and  Orlov  [20]. 


In  Chapter  III  we  derive  two  asymptotic  eq'  ations,  valid  in 
different  regions,  which  describe  the  behavior  of  Korteweg-deVries 
solitary  waves  in  slowly  varying  media.  In  the  core  we  use  quasi- 
stationarity  to  obtain  the  asymptotic  solution.  However,  in  the 
leading  tail  quasi-stationarity  is  not  appropriate,  thus  we  use  the  more 


standard  multiple  scale  approach. 

We  next  show,  in  Chapter  III,  that  the  problem  of  a  KdV  solitary 
wave  propagating  through  slowly  varying  media  is  similar  in  some 
respects  to  the  accelerating  piston  problem  in  gas  dynamics.  The  method 
of  characteristics  shows  that,  for  certain  physically  interesting 
problems,  characteristics  cross  in  the  leading  tail  causing  a  caustic  to 
form  which  predicts  a  region  of  triple  valuedness.  By  examining  the 
local  behavior  near  the  point  of  first  crossing,  we  find  it  useful  to 
introduce  a  similarity  variable.  The  breakdown  of  the  asymptotic 
expansion  leads  to  rescaling.  The  new  scales  show  that  the  leading 
order  solution  near  the  first  crossing  of  characteristics  satisfies  the 
diffusion  equation,  even  though  the  wave  could  be  purely  dispersive 
without  dissipation.  Using  Laplace's  method  for  the  asymptotics  of 
exponential  integrals,  we  show  that  the  solution  does  not  become  triple 
valued  (break),  but  instead  rapid  transitions  occur  from  one  exponen¬ 
tially  decaying  tail  to  another,  a  wave  number  shock. 

In  Chapter  III  we  show  that  a  penumbral  caustic  is  generated  when 
the  first  crossing  of  characteristics  occurs  along  the  characteristic 
eminating  from  the  initial  core  position.  Then,  in  Chapter  IV,  we  show 
for  the  other  two  cases  (first  crossing  along  a  characteristic  eminating 
from  the  moving  core  at  a  later  position  or  eminating  from  the  initial 
leading  tail)  that  a  cusped  caustic  is  generated.  Irt  all  three  of  these 
cases  we  show  that  the  predicted  triple  valuedness  actually  corresponds 
to  the  propagation  of  a  wave  number  shock  in  the  leading  tail. 
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CHAPTER  II 


■%*  -„' 
i  "  *> 


DERIVATION  OF  THE  EQUATIONS 

2.1  Introduction 

In  this  chapter  we  seek  partial  differential  equations  (pde's), 
using  perturbation  techniques,  which  describe  the  behavior  of  slowly 
varying  solitary  waves  and  solitons  [1]  moving  in  inhomogeneous  media. 

Since  the  pioneering  work  of  Gardner,  Greene,  Kruskal  and  Miura 
[4]  in  1967,  it  has  become  well  known  that  the  constant  coefficient 
Korteweg-deVries  (KdV)  equation  has  a  family  of  exact  N-soliton 
solutions.  These  solutions  can  be  obtained  using  the  inverse  scattering 
transform  technique  [1,29]  which  was  derived  for  exactly  integrable 
equations.  The  variable  coefficient  KdV  equation,  which  in  general  is 
not  exactly  integrable,  describes  many  interesting  physical  phenomena 
[6,  11,  12,  17,  21,  25,  26]  and  has  been  discussed  extensively  in  the 
literature.  Grimshaw  [5]  and  Johnson  [6]  consider  the  case  of  a 
solitary  wave  propagating  on  the  surface  of  water  of  variable  depth.  We 
make  a  minor  generalization  to  their  work  by  allowing  the  variable 
coefficients  to  depend  on  space  as  well  as  time 

Using  a  multiple  scale  approach,  Grimshaw  [5]  showed  that  the 
asymptotic  expansion  breaks  down  as  one  gets  too  far  from  the  peak  or 


k.v. 


K\ 

V.'w 


m 


center  of  the  wave.  We  thus  consider  two  cases,  first  the  "core"  or 
thin  region  near  the  center,  and  then  the  "leading  tail"  or  semi- 
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(2.5) 


-wU  +  vkUU  +  Ak  U 

n  n  nnn 


-e(UT  +  vUUv  +  3AkkvU  +  3AkvllUv  +  yU) 

T  X  X  nn  X  Xnn 

-  e2(\kxxUn  +  3XkxUXn  +  3AkUXXn)  -  e3(AUXX)()  . 

To  solve  this  equation,  we  expand  U,  k  and  w  in  powers  of  e 
U  =  Ug  +  eU^  +  e“U 2  +  ...  , 

2 

k  =  kQ  +  ek1  +  e  k2  +  ...  , 


(2.6) 


w  =  Wg  +  ew^  +  e  w^  +  ... 

Thus,  to  leading  order,  equation  (2.6)  becomes 

■w0U0n  +  vk0U0U0n  +  Xk0U0nnn  =  0  *  (2,7' 

The  well  known  solitary  wave  solution  for  equation  (2.7)  is 

UQ  =  A  sech28n ,  (2.8) 


v  k  ’ 
0 


8  =  (-V/2  - 

4Ak 

K0 

which  agrees  with  Grimshaw  [5].  Continuing  to  higher  orders  shows  that 
the  expansion  for  U  breaks  down  when  |n |  >  0(1).  We  thus  define  the 
core  by  n  =0(1)  (assuming  g  =  0(1)),  since  otherwise  Ug  is  exponen¬ 
tially  small.  The  center  of  the  core  is  n  =  0.  We  take  the  position  of 
this  center  to  be  an  unknown  function  of  time 


X  =  X  (T)  =  X  (T)  +  EX  (T). 
c  0  C1 


(2.9) 


Thus,  at  any  given  time,  we  find  from  equation  (2.4. a)  that  the  fast 
phase  expanded  around  its  center  is  given  by 
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n  =  kQ(XcJ)[X-Xc]/e, 


(2.10) 


to  leading  order.  Thus  the  core  occurs  where  X-Xc  =  0(e)  (where  (2.10) 

is  valid),  or  to  within  order  £  of  the  path  X  =  X  (T).  Now  using  "  A" 

C0 

to  denote  evaluation  at  X  =  X  ,  we  find  that  the  leading  order  solution 

C0 

(2.8)  in  the  core,  is  given  by 


it  -  3w°  h2rfW°  i1/2  X"Xci 

Un  =  — t-  sech  L(-rx)  - J  • 

Ok0  4Xk0 


(2.11) 


This  is  a  quasi-stationary  solution  [18].  We  redefine  the  fast  phase 


to  be 


From  (2.4)  this  choice  of  phase  yields 
k  =  1  , 


w  =  X '  (  T )  , 
c 


(2.12) 


(2. 13. a) 


(2.1 3. b) 


automatically  satisfying  (2.5).  We  will  shortly  determine  the  unknown 


w(T).  In  terms  of  it  the  center  of  the  core  is 
Xc(T)  =  Xc(0)  +  fJQ w(T)dT  , 


(2.13.c) 


where  Xc( 0 )  is  the  initial  center.  It  is  thus  easier  to  directly  solve 
(2.1)  in  the  core  using  a  quasi-stationary  approach  with  two  scales, 
rather  than  the  standard  slowly  varying  approach  with  three  scales: 


u(x,t)  =  V(n ,T) • 


We  expand  the  coefficients  around  X  =  X  to  obtain 

a  ~  c0 

-wV  +  vW  +  XV  = 


-e[VT  +  vY(n  +  X„  )vvw  +  XY(n  +  X„  )V  +  YV] 
t  X  c1  n  X  c1  nnn 


+  0(e )  . 


(2.14) 


(2.15) 


Thus  the  leading  order  equation  is 

-wnVn  +  vVnVn  +XVn  =  0  , 
0  On  0  On  Onnn 


(2.16) 


with  the  solution 


VQ  =  ^sech2[(^L)1/2 


X  -  /gV(T)dT 


(2.17) 


Here  we  have  taken  the  center  of  the  core  to  be  initially  at  X  =  0.  Note 

A  .  . 

that  the  sign  of  w^  equals  the  sign  of  X  and  hence  equations  v2.13.by 
and  (2.17)  show  that  the  core  must  move  in  the  direction  of  the  sign  of 

X.  Continuing  to  the  next  order  we  obtain 

-w nV1  +  C(VnVj  +  XV.  =  w  Vn 
0  in  o  1  n  1 nnn  i  On 


-  (V0T  +  GX(ri  +  XCl)VoV  ^  +  XC;L)V0nnn  +  'V  • 


(2.18) 


We  now  define  the  linear  operator,  L,  and  its  adjoint,  L  ,  as  did 
Kodama  and  Ablowitz  [18]: 

Ly  i  -V„  *  v(v0y)n  *  1y  . 


L  y  =  W  y  -  vVny  -  Xy  , 
On  On  nnn 


so  that 

ov-h  -  viL\;d"  = 

[+x(Vlnn  -  Win  *  VmV  +  *V0V1  '  Wl]l  • 

Thus  and  its  first  two  partial  derivatives  with  respect  to  n  will  be 
well  behaved  if  the  following  capatibility  condition  is  satisfied 

Ov ovoT  *  :x(n+xc1,v§vo,1  ♦  *vo!d''  --  0  •  (2-19! 

Using  equation  (2.17)  to  evaluate  (2.19)  we  obtain 


‘1/2w3/2 


dT  v  "2 

v 


O  xV2w3/2 

X  .96  X  ..  „:w  W0  \ 

~  (  5  *  w0  "  2Y/(  ~ 1  ^  ’ 

V 


(2.20) 


sJl  JL.  .  “V  A*.  »/■  *■  *-*».-»■  \m.r.  «  •*>  yj\  V  .  t.A,  -•’> 
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which  can  be  solved  for  Wg(T)  from  an  initial  condition  w(0)  (we  assume 

w(Q)  independent  of  £  so  that  Wg(Q)  =  w(0)).  Thus,  the  leading  order  core 

solution,  (2.17),  is  completely  determined  using  (2.13.b). 

Equation  (2.20)  appears  to  be  a  first  order  ordinary  differential 

equation  (ode)  with  variable  coefficients.  However,  in  general  the 

position  of  the  core  depends  on  the  unknown  w^  via  equation  (2.13.c). 

We  see  that  w_  =  X'  and  wi  =  X"  .  Thus  (2.20)  is  actually  a 
u  c0  u  cQ 

second  order  ode  with  variable  coefficients.  To  solve  it  we  must  know 

the  initial  conditions  X  (0)  and  X'(0).  If  v  and  X  are  functions  of  T 

c  c 

only,  then  (2.20)  has  the  solution 

,tN  ,nxr  X(0)v4(T),1/3  r  4.T  ...  .T, 
wn(T)  =  w(0)[  -r - ]  exp[-T/nY(T)dTJ  , 

U  v4(0)X(T)  5  0 

as  shown  by  Ablowitz  and  Segur  [1], 

Vie  note  here  that  if  the  perturbation  analysis  using  equation 
(2.3)  is  carried  forward  without  assuming  quasi-stationarity,  we  can 
obtain  a  hyperbolic  system  of  pde’s,  in  X  and  T,  for  the  wave  number  and 
amplitude,  as  did  Grimshaw  [5].  We  must  be  careful  to  realize  that 
these  pde's  are  only  valid  in  the  core;  X  is  in  a  very  thin  region  near 
the  center.  From  equation  (2.10)  we  see  that 
X  s  Xc(T)  +  0(e). 

In  this  thin  region,  X  is,  to  leading  order,  a  function  of  T,  thus 
reducing  the  pde's  to  ode's.  If  we  fail  to  show  proper  care,  and  use 
the  method  of  characteristics  to  solve  the  hyperbolic  system,  after  much 
tedious  algebra  and  some  complicated  matching,  we  find  that  the  solution 
fails  to  focus  in  the  core  region.  This  is  the  same  result  we  have 
obtained  using  quasi-stationarity! 


*>  w.  *-*  “-r  ->  V  *.v  >  */-.* -V-V-  J  nm^  V-  *  L  -i-  «  - VL  *  -4  K*'- 
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2.3  The  Leading  Tail 

The  semi- inf inite  regions  on  either  side  of  the  core  we  refer  to 
as  tails.  The  tail  on  the  side  of  the  core  in  the  direction  of  its 
motion  we  call  the  leading  tail.  We  call  the  other  one  the  trailing 
tail.  There  is  much  discussion  in  the  literature  concerning  the 
trailing  tail  of  slowly  varying  solitary  waves  and  the  development  of  a 
"shelf-like"  structure  [13,  16,  17,  18,  19,  24,  25,  28].  There  appears 
to  be  very  little  discussion  concerning  the  crossing  of  characteristics, 
which  leads  to  the  formation  of  a  caustic,  for  the  leading  tail  [5].  It 
is  this  problem  we  are  concerned  with  here. 

In  the  leading  tail  the  quasi-stationarity  assumption  is  not 
valid,  since  the  tail  is  not  localized  in  space.  However,  due  to  the 
decay  as  of  the  solution,  (2.17),  the  nonlinear  term  in  equation 

(2.1)  is  exponentially  small.  Thus,  as  is  justified  by  the  method  of 
matched  asymptotic  expansions,  we  seek  a  slowly  varying  solution  to 

u.  +  Xu  =  - eyu  .  (2.21) 

t  xxx 


We  expect  the  leading  tail  of  the  solitary  wave  to  be  exponen¬ 
tially  decaying.  A  steady  traveling  decaying  wave  solution  would  be 

._KX-flT 


similar  in  mathematical  form  to  an  oscillatory  dispersive  wave.  However, 
(2.21)  has  slow  variations  including  a  slow  decay  or  growth  introduced 
by  the  term  with  .  Thus  we  use  the  multiple  scale  approach  and  intro¬ 
duce  terminology  similar  to  that  used  fir  oscillatory  dispersive  waves: 
the  fast  phase,  H,  slowly  varying  wave  number,  K,  and  slowly  varying 
frequency,  .  They  satisfy  the  conditions 


f 

1 •  *  *  r 


^  * 
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Hx  =  e(X,T)/e, 


(2.22) 


Hx  =  ex  H  K  , 


(2. 23. a) 


=  8r  -  , 


( 2 . 2  3 .  b ) 


so  that  conservation  of  waves 


k,  +  nx  =  o  , 


also  holds  in  the  leading  tail.  Using  the  chain  rule  on 
u(x,t)  =  U(X,T,H), 
equation  (2.21)  becomes 

-ftUH  +  '^X?UHHH  =  -e(Ut  +  3XKKxUhh  +  3AK2UXHH  +  YU) 


(2.24) 


-e2(AKXxUH  +  3XKXUXH  +  3aKUxxh)  -  e?(xUm)  . 


(2.25) 


The  leading  order  equation 


-VOH  *  1K0U0HHH  =  0  - 


(2.26) 


has  the  decaying  solution  (as  H  assuming  A>0  or  a  right  going 


Un  =  A(X,T)exp[-(-~)1/2H  ]  . 

K 

The  next  order  equation  is 


(2.27) 


-VlH  *  iK0U1HHH  =  !,,U0H-»K0I<1U0HHH  ’  U0I 


“3XKokoxuohh  ■  3XKouoxhh  “  yUo  * 


(2.28) 


The  right  hand  side  of  this  equation  produces  secular  terms  for 

(terms  which,  for  large  H,  cause  the  asymptotic  expansion  to  break  down) 

unless 


Q0  =  XK0  ’ 


(2.29) 


-  /\-  v 


V  V  V  '  .*  V  V  ’ .*  *.>  V  *.'•  '.k  V  V  V  W  *  •  V.  '  ■  '  ■  -  . •  **  -  *  .  ■  •  *■  .  “  *  .  •  .  "  .  ■  _ .  *  •  '  .  ■  •  ■  -  ' 
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At  +  3AKqAx  =  -( y  +  3XKqKqx  -  3XKqK1  +  fl^A. 

Using  (2.29)  and  (2.30)  in  (2.28)  we  can  solve  to  obtain 


U1  =  B(X , T)exp[-H]  . 


(2.30) 


(2.31) 


From  equations  (2.27)  and  (2.29),  and  (2.31)  we  see  that  it  is  not 
necessary  to  expand  H,  the  phase,  since  the  order  e  and  higher  order 
terms  can  be  absorbed  into  the  amplitudes.  Thus  equations  (2.29)  and 


(2.30)  become 

0  =  XK3  , 


Ay  +  3XKZAX  =  -(y  +  3XKKx)A. 


(2.32) 


(2.33) 


Equation  (2.32)  is  the  "dispersion  relation"  for  exponential  tails  for 
the  linearized  KdV  equation.  Using  this  relation,  equation  (2.24)  for 


the  conservation  of  waves  becomes 
Kt  +  3XK2Kx  =  -XXK3  , 


<"2.34) 


which  must  be  solved  along  with  equation  (2.33).  The  leading  tail 
solution  must  then  be  matched  to  the  core  solution.  Under  physically 
reasonable  conditions,  we  will  show  in  Chapter  III,  using  the  method  of 
characteristics  for  equation  (2.34),  that  characteristics  will  cross 
causing  a  caustic  to  form.  Then  we  will  consider  this  caustic  formation 
and  its  interpretation. 


CHAPTER  III 


ACCELERATING  PISTON  PROBLEM 

3,1  Introduction 

In  their  now  classic  monograph  [3],  Courant  and  Friedrichs 
briefly  discuss  a  gas  dynamics  problem  in  which  a  compression  wave  is 
produced  by  a  uniformly  accelerated  piston.  Figure  3.1  shows  the 
envelope  formed  by  the  crossing  of  straight  characteristics  which  de¬ 
scribe  this  problem.  In  this  chapter  we  show  that  the  moving  core  of  a 
slowly  varying  solitary  wave  acts  like  the  piston  of  this  gas  dynamics 
problem.  In  Section  3.6  we  show  that  the  constant  coefficient  problem, 
equation  (2.1)  with  constant  coefficients,  has  straight  characteristics 
which  for  an  accelerating  wave  form  the  same  envelope  as  the  uniformly 
accelerated  piston  (Figure  3.1).  If  on  the  other  hand  the  wave  is 
decelerating,  then  a  rarefaction  wave  is  produced  as  with  a  uniformly 
decelerated  piston  (Figure  3.2).  In  Section  3.2  we  show  for  the  more 
general  case  of  variable  coefficients  that  an  envelope  may  or  may  not 
form  depending  on  the  medium  and  on  the  motion  of  the  wave.  The  usual 
case  of  an  accelerating  wave  with  variable  coefficients  is  illustrated 
in  Figure  3.3,  where  the  characteristics  are  now  curved  due  to  the 
variable  coefficients. 

We  also  show  in  Section  3.2  how  the  characteristics  arise  and 
can  lead  to  a  focusing  solution.  Then,  in  Section  3.3,  we  utilize  a 


Figure  3.1  Envelope  of  straight  characteristics 
for  a  compression  wave  produced  by  a 
uniformly  accelerates  piston  [3], 


Figure  3.2  Straight  characteristics  of  the  rarefaction 
wave  produced  by  a  uniformly  decelerated  piston  [3]. 


Figure  3.3  Envelope  of  the  curved  characteristics  for 
the  leading  tail  of  a  slowly  varying  solitary  wave 
propagating  on  the  surface  of  water  of  variable  depth. 


similarity  variable  to  determine  equations  for  the  envelope  ^boundary  of 
the  region  of  multiple  characteristics),  and  to  obtain  a  "fundamental 
folding  equation".  Examining  the  local  behavior  near  first  focusing 
leads  us  to  new  scales.  We  show,  in  Section  3.4,  how  these  new  scales 
relate  the  solution  to  a  solution  of  the  diffusion  equation.  In  Section 
3.5,  we  use  Laplace's  method  for  the  asymptotics  of  exponential 
integrals  to  determine  the  asymptotic  behavior  of  our  solution  near 
first  focusing.  We  find  the  solution  to  be  a  single  exponential  tail 
everywhere  except  in  a  relatively  thin  neighborhood  of  a  "shock  path", 
where  the  solution  makes  a  transition  from  one  exponential  tail  to 
another.  In  this  neighborhood  the  solution  is  the  sum  of  two  tails.  The 
"shock  path"  is  actually  the  location  in  space  time  of  a  rapid  change  in 
wave  number  or  a  "wave  number  shock".  This  idea  of  a  wave  number  shock 
was  first  discussed  for  slowly  varying  solutions  of  reaction-diffusion 
equations  by  Howard  and  Kopell  [10].  Finally,  in  Section  3.6,  we  work 
through  an  example  of  our  analysis  using  constant  coefficients. 


3.2  A  Focusing  Solution 


In  the  last  chapter  we  determined  the  leading  order  solution  for 


both  the  core 


V0  =  3 


—  sech2 


[(¥)'"]  ■ 


and  the  leading  tail  (right  go:  ng  with  X  >  0) 
Uq  -  A  exp[-H]  . 


0.1 ) 


(3.2) 


Using  (2.1 3.b)  and  the  given  coefficients  v  ,  X  and  y ,  we  can  solve 
equation  (2.20)  to  completely  determine  Vg,  and  thus  to  leading  order 


the  path  of  the  core,  X  (T).  We  now  seek  the  behavior  of  IL  by  solving 

c0  U 
the  system  of  pde's  given  by  equations  (2.33)  and  (2.34).  First  we 

match  tig  and  Vg  using  the  basic  matching  principle  (see  [15,20]).  We 

write  Vg  in  terms  of  X  and  T,  Ug  in  terms  of  n  and  T,  and  expand  for 

£  -*■  0.  Thus  we  have 


x '  X' 

cn  c_  1/9  X-X„ 

V.  'V  12  —  exp[-(— ;1/2  - £]  , 

0  c  5 


Un  -v  A  exp[-K  - -]  , 

u  z 

where  "  * "  again  means  evaluation  at  X  =  X  .  Now  the  matching  of  Ug  and 

co 


Vg  yeilds  the  matching  conditions 


A  =  12 


(3.3) 


■(-r) 


Co  \  1/2 


(3.4) 


Next  we  use  the  method  of  characteristics  to  solve  the  leading  tail 
equations  (2.33)  and  (2.34): 

At  +  3AK2Ax  =  (y  +  3AKKX)A  , 


Ky  +  3 AK2Kx  =  -AXK3  . 


T  aking 


UA  _  X'  v1- 
■pr  -  AaK  , 


(3.5) 


.  *.  %  *.  *.  \  *.  *-  *.  *.  *.  •.  •.  •  %  -v  *.  *.  A  %  \  .%  .*•  A.V 
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we  see,  using  equation  (3.4) ,  that  near  the  core  the  characteristic 
velocity  for  the  leading  tail  is  three  times  the  core  velocity.  This 
means  that  disturbances,  in  the  leading  tail  near  the  core,  will  move 
away  from  the  core  into  the  tail.  Thus,  the  moving  core  acts  as  a 
space-like  boundary  for  the  leading  tail.  The  problem  of  solving  the 
leading  tail  equations  subject  to  a  specified  core  (amplitude  and  phase 
is  well-posed;  there  is  a  unique  leading  tail  matching  to  the  core.  On 
the  characteristics  our  system  of  pde's  reduces  to  a  system  of  ode's 

HA 

■gy  =  -(y  +  3xKK^)A  ,  (3.6) 

and 


£  ■  V 


(3.7) 


Equations  (3.3)  and  (3.7)  imply  that  the  characteristics  are  given  by 
X  =  F(T;£)  ,  (3.8) 

where  £  is  the  parameterization  of  time  along  the  path  of  the  leading 
order  core  (at  T  =  £,  X  =  X  (£)): 


0 

F(c ;c)  =  X  (£),  £  >_  0  . 

0 


(3.9) 


We  call  the  crossing  of  characteristics  focusing,  and  the  locus  of  all 
such  focusing  points  is  a  caustic.  The  caustic  is  obtained  by  simul¬ 
taneously  solving  (3.8)  and 

F  (Tf;£)  =  0  .  (3.10) 


The  times,  T  =  T^,  for  which  focusing  occurs,  due  to  the  moving  core, 
are  obtained  by  solving  equ&r.ion  (3.10)  for  T^: 

Tf  =  Tf(£)  .  (3.11) 


Focusing  occurs  ahead  of  the  core  if 


Tf(£)  >  £. 


The  corresponding  positions  at  which  focusing  occurs  are  then  determined 
from  equation  (3.8): 

Xf  =  F(Tf;5)  .  (3.12) 

Equations  (3.11)  and  (3.12)  are  the  parametric  representation  of  the 
caustic  due  to  the  moving  core. 

Thus  far  we  have  considered  only  those  characteristics  eminating 
from  the  path  of  the  core  (see  Figure  3.4).  We  now  consider  the  charac¬ 
teristics  due  to  initial  conditions  which  eminate  from  the  positive 
X-axis  (see  Figure  3.5).  Equations  (3.5)  and  (3.7)  imply  that  these 
characteristics  are  given  by 

X  =  G(T;;)  ,  (3.13) 

where  S  is  the  parametrization  of  initial  position  (at  T=0,  X=0: 

G(0*,C)  =  ?  ,  ?  >  0. 

Analyzing  these  characteristics  as  above  we  see  that  they  too  can  focus. 
We  do  not  pursue  the  case  where  neither  (3.8)  nor  (3.13)  leads  to 
focusing  since  there  is  then  no  difficulty  with  our  asymptotic 
expansions.  We  proceed  under  the  assumption  that  there  is  focusing  due 
to  the  movement  of  the  core  only,  and  will  consider  the  other 
interesting  case  of  focusing  due  to  initial  conditions  in  Chapter  IV 
(see  also  [8]). 

Returning  to  the  characteristics,  (3.8),  we  determine  the  time  of 
first  focusing,  T=T^  ,  by  minimizing  (3.11)  for  £>0.  There  are  two 
possibilities.  First  we  can  have 

T1  =  Tf (0) ,  (3.14) 

with 


T1C0)  >  0  . 


(3.15) 


**  / 
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In  this  case  first  focusing  occurs  along  the  characteristic,  5=0, 
eminating  from  the  origin;  focusing  due  to  later  movement  of  the  core 
happens  at  a  later  time  (thus  giving  rise  to  a  penumbral  caustic,  see 
Fig.  3.4  and  [20]).  The  other  possibility  is 

T1  =  Tf(51)  ,  (3.16) 

with 

T^)  =  0  (usually  with  T"(^)  >  0)  .  (3.17) 

In  this  case  first  focusing  occurs  along  the  characteristic, 5  =  5^  >  0; 
focusing  due  to  earlier  and  later  movement  of  the  core  occurs  later 
(thus  giving  rise  to  a  cusped  caustic).  We  consider  the  first 
possibility,  equations  (3.14)  and  (3.15),  in  the  remainder  of  this 
chapter.  The  second  possibility,  equations  (3.16)  and  (3.17),  will  be 
covered  in  Chapter  IV  along  with  focusing  due  to  initial  conditions. 


3.3  First  Focusing  from  the  Initial  Core  Position  (Local  Behavior) 


We  now  consider  the  case  where  first  focusing  occurs  due  to  the 


immediate  movement  of  the  core  from  its  initial  position,  that  is 
equations  (3.14)  and  (3.15)  apply  (of  course  the  time  to  first  focus  is 


finite,  T^(0)>0).  We  approximate  the  equation  for  the  characteristics 
near  first  focusing  by  Taylor  expanding  (3.8)  around  5=0  and  T=T^  =  Tjr(0). 
In  this  way  we  obtain  the  "fundamental  folding  equation": 


X-X1-Ft(T1;0)(T-T1)  -  ^Fu(T1;0)(T-T1)2  = 


FTc(T1;0)(T-T1)5  +  If  ^  ;0)  52 , 


(3.18) 


a  cut  off  quadratic  in  5.  From  equations  (3.10)  and  (3.11)  we  see  that 
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dc  r5(rf=5;  *  Fu(rr5:Tr( e:  -  =  0  •  !3-,9) 

Thus  if  TJ.(0)  =  0,  then  F  (T^O)  =  0  and  we  would  need  to  include 
cubic  terms  in  the  approximation  (3.18).  The  analysis  for  the  cubic 
will  be  done  in  Chapter  IV.  We  now  proceed  under  the  assumption  that 
T^(0)  >  0.  Equation  (3.18)  is  valid  for  e,  >_  0,  since  the  parame¬ 
terization  for  the  position  of  the  core  begins  at  £  =  0.  We  notice  from 
equation  (3.18)  that  the  change  of  variable 


£  =  (T-T1)f(S), 


(3.20) 


shows  that  f(S)  is  a  function  of  the  similarity  variable  S,  where 


S  i  [X-X1  -  C  ( T-T1 )  -  C(T-T1)2]/(T-T1)2, 

with  the  constants  C  and  C  defined  by 

0 

V  FT(T1;0)  ’ 


(3.21) 


C  =  -ifTT(Ti;0)  . 


Thus  the  approximation  for  the  characteristics,  (3.18),  becomes 
af^S)  +  bf(S)  =  5  , 

where  the  constants  a  and  b  are  defined  by 
a  =  j  F  (TlS0), 


(3.22) 


b  £  FT€(T1 ;0). 


From  equation  (3.19)  we  then  have 


a  =  -1  T^(0)b. 


(3.23) 


Now,  for  T<T,  ,  X-X,  <  C  (T-TJ  +  C(T-T,r,  and  5  >  0,  we  see  from  (3.21) 
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that  3  <  0.  From  (3.22)  we  have 

f(S)  =  [-b  ±  Jq-  +  4aS]/2a  ,  (3.24) 

and  (3. 20)  implies  f(S)  <  0.  In  this  region  of  space  time  there  is  no 
crossing  of  characteristics,  so  there  must  be  a  single  characteristic 
for  each  S.  From  (3.24)  we  see  that  in  order  for  f(S)  to  have  a  single 
negative  value,  and  thus  from  (3.20)  for  £  to  have  a  single  positive 
value,  we  must  have  aS  >  0  or  a  <  0.  Equation  (3.23)  then  yields  b  > 
0.  Examining  equations  (3.20)  and  (3.24)  for  T  >  shows  that  these 
signs  for  a  and  b  yield  the  required  two  positive  values  of  f(S)  in  the 
appropriate  region  of  crossing  characteristics.  Physically,  b  >  0  means 
the  velocity  of  the  leading  tail  is  increasing,  with  respect  to  the 
position  of  the  core,  at  first  focusing  (£  =  0  and  T  =  T. ) . 

Since  is  the  time  of  first  focusing,  all  focusing  occurs 

for  times  T>T^  ,  and  thus  the  envelope  of  crossing  characteristics 
is,  given  by  equation  (3.24),  the  boundary  of  the  region  where  f(S) 
has  two  positive  values: 


S  =  0,  T  >  T1  , 

and 

h2 

S  =  -  r-  ,  T  >  T.  . 

4a  ’  —  1 

The  neighborhood  of  first  focusing  can  be  divided  into  three  regions 
according  to  the  number  of  characteristics  passing  through  each  point 
(see  Figures  3.3  and  3.6).  Region  I  is  given  by 


X  <  X1  +  C  (T-T1)  «•  Cd-T^' 


and  has  a  single  characteristic,  eminating  from  the  moving  core,  passing 
through  each  point,  i.e.  f(S)  has  one  value  satisfying  both  (3.20)  and 
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(3.22).  Region  II  is  given  by 

*1  *  Cg<T-V  "  C(T-V2  <  X  <  X,  *  cg<  r-T, )  ♦  (c  -  f>V 

This  region  has  two  characteristics,  eminating  from  the  moving  core, 
passing  through  each  point,  i.e.  f(S)  has  two  values  satisfying  both 
(3.20)  and  (3.22),  and  a  third  characteristic,  eminating  from  the 
initial  conditions,  also  passes  through  each  point.  Region  III  is  given 


X  >  X1  +  C  (T-T1)  +  C C T— T 1 )  ,  T  <  T1  , 


X  >  X1  +  C  (T-T^  +  (C  -  -55)(T-T1)2,  T  >  Tr 

This  region  has  a  single  characteristic,  eminating  from  the  initial 
conditions,  passing  through  each  point  and  no  value  of  f(S)  satisfying 
both  (3.20)  and  (3.22).  Figure  3.6  shows  these  regions,  and  Figure  3.7 
gives  the  evolution  of  the  number  of  characteristic  values  at  each  X  due 
to  the  moving  core. 

We  now  seek  the  behavior  of  K  and  A  near  first  focusing,  X=X^  and 
T=T^.  From  equation  (3.8)  we  have 

K ( X , T)  =  K[F(T;C),T]  =  K(T;0  , 
and  expanding  in  a  Taylor  series  around  T=T^  and  5=0 


K  =  R(T1  ;0)  Rj(T1  ;0)(T-T1 )  +  K  (1^0)5  +  ...  . 


Substituting  (3.20)  into  this  expression  yields 

K  =  K ( T 1 ; 0 )  +  [Kt(Ti;0)  +  K . ' T 1  ;  0 )  f ( 5 ) ] ( T - T  ^ )  +  ...,  (3.25) 


and  thus  as  T  -*  T^  (keeping  S  fixed  so  that  5  -*•  0) 


Kx  ^  R  (T, ;0)f,(S)/'T-T1  )  . 


(3.26) 


From  equations  (3.5)  and  (3.18)  we  have 


3XK2  'v  C  +  b^  , 

g 

and  thus 

3(XXK2  +  2XKKX)  'v  bCx  .  (3.27) 

For  well  behaved  X,  (3.25),  (3.26)  and  (3.27)  yield 
6XK«X  %  b5x  , 
and  from  (3.18) 


Thus, 


6XKKy  'v  't  i'  - 
X  bvT-T^ ;  +  2a£ 

equations  (3.6)  and  (3.28)  yield 


(3.28) 


dA  1  bA 

dT  ^  "  2  dTT-T  )  '+  2ar 

or 

A  *  S( C) |b(T-T1)  +  2aci‘1/2  , 

where  A(5)  is  determined  by  the  matching  condition  (3.3).  Using 
(3.20)  we  now  have 

A  *  5( r)|[b  +  2af(5)](T-T1) |"1/Z  .  (3.29) 

The  singular  behavior  of  A  at  first  focusing  suggests  the  use  of 
different  scales,  which  we  introduce  in  the  following  section. 


3.4  A  Relation  to  the  Diffusion  Equation 
From  equation  (2.25)  we  obtained  the  leading  order  equation  (2.26 
with  solution 


Un  =  A( X  , T )  expC-H]  . 


'5.30) 


Near  first  focusing  A  is  singular  as  given  by  (3.29).  We  will  show 
higher  order  terms  in  the  expansion  are  more  singular.  At  the  next 
order  we  obtained  the  solution 


U1  =  B(X,T)  exp[-H]  .  (3.31) 

The  next  order  equation  is  very  important: 

-  ftU2H  +  AK  U2HHH  =  -yU1  U1T  xKxxU0H 

3XKKXU1HH  -  3AKxUqxh  -  3XK  XHH  -  3aKUQX)(H  (3.32) 

We  see  that  the  right  hand  side  produces  secular  terms  for  U2  unless 

Bt  +  3AK2Bx  =  -(y  -t-  3AKKX)B  +  AK X)(A 

+  3AK)(Ax  +  3AKAxx  .  (3.33) 

The  asymptotic  behavior  of  A,  (3.29),  leads  us  to  assume 

B  ^  1  T-T1  |P  T(S) .  (3.34) 

Substituting  (3.23),  (3.29)  and  (3.34)  into  (3.33)  yields  P  =  -  7/2. 
Continuing  this  process  for  higher  orders  shows  that  the  expansion  for  IJ 
breaks  down  when 

T-T1  =  0(e1/3)  , 

which  suggests  the  scaling 

t  =  e"1/3(T-T1 )  =s2/3(t-t1)  .  (3.35) 

In  order  to  keep  the  similarity  variable,  S,  order  one,  we  must  take 

Z  =  c"2/3[X-X1-Cg(T-T1)-C(T-T1)2] 

=  e  1/3[x-x.-C  (t-t'-eCU-t.)2]  .  (3.36) 


Equations  (3.35)  and  (3.36)  are  the  appropriate  independent  variables  in 
a  neighborhood  near  first  focusing.  We  note  that  the  scales  of  these 
variables  are  intermediate  to  the  other  two  scalings,  large  compared  to 
the  scales  of  the  solitary  wave  and  small  compared  to  the  scales  of  the 
slowly  varying  media.  Using  equations  (2.22),  (2.23)  and  (2.32),  near 
first  focusing  we  have 

H  =  ■1[0(X1  ,T1  )  +  K(X1,T1)(X-X1)  -  x(X1  .T^K3^  J^U-T,,)  +  ...] 

=  1  [0(1)-K(1)X1  +  X(1)K(1)3T.]  +  K(1)[x-A(l)K(1)2t]  +  ...  , 
e  1  1  (3.37) 

where  the  superscript  "(1)"  means  evaluated  at  X=X^  and  T=T^ .  Since, 

from  (3.29),  the  amplitude  is  singular  at  first  focusing,  equations 

(3.30)  and  (3.37)  motivate  us,  using  the  intermediate  scales  (3.35)  and 

(3.36),  to  seek  a  solution  of  the  form 

u  =  4>(Z,T)exp[-K(1'(x-\(1)K(1)  t)]  ,  (3.38) 


where  the  large  constant  in  (3.37)  will  be  part  of  <j>(Z,i).  This  is  a 
different  multiple  scale  approach  than  the  one  we  began  with  in  Section 
2.2  (equation  (2.3)).  Thus  we  have 


u ( x , t )  =  u(x,t,Z,r)  , 


(3.39) 


ut  =  [X(1)K(1)  $  -  e1/3Cg4>z  +  e2/\x  -  2e4/3C( t-^  )#z] 


cp[-K(1)(x-X(1)K(1)  t)]  , 


ii  r  a  ^  x  1  /3^(1 ) 2  ,  2/3w, ( 1 )  ^  A  -i 

=  L-K  +  3e  K  -  3e  K  4>zz  +  e<J>zzzJ 


f  1  1  M'  (  1  \  2 

<p[-Kv1/(x-Av1/Kn'  t)  ]  . 


Substituting  these  values  into  equation  (2.21)  we  have 

,.(1V!1).  1/3  re-,  (1  )  /  \  ,(1),.  ,  .  1„(  1  )3 

=  3  a  K  +  e  tUx  (x-x^)  +  At  ( t—  t ^ )  ]K  p 


-  -\(1)1>zzz  -  Y(1l>  +  H.O.T. 

Near  first  focusing  the  leading  order  term  for  satisfies  the  diffusion 
equation 


*o  =  5X<1>k(1Vz  • 


(3.40) 


We  showed  in  Chapter  II  that  for  the  right  going  wave,  which  we  have 
considered  here,  X  >  0,  and  hence  (3.40)  is  the  well-posed  diffusion 
equation.  It  is  extremely  interesting  to  note  that  even  if  y  =  0  we 
still  easily  obtain  the  diffusion  equation  from  the  KdV  equation  even 
though  the  KdV  equation  is  dispersive  with  no  dissipation!  If  we 
momentarily  consider  the  crossing  of  characteristics  for  the  trailing 
tail  of  a  right  going  wave,  then  the  analysis  is  nearly  identical  to  the 
above.  However,  equation  (3.38)  becomes 


u  =  ^  exp[K(1)(x-X(1)K(l)  t)], 

and  leads  to  the  backwards  diffusion  equation,  which  is  ill-posed.  In 
this  case  the  solution  is  grossly  unstable,  and  therefore  the  assump¬ 
tion  of  a  slowly  varying  exponentially  decaying  tail  is  not  valid  for 
the  trailing  tail  of  a  right  going  wave.  In  fact,  the  trailing  tail  is 
known  to  not  be  exponentially  decaying.  Ko  and  Kuehl  [16,  17],  Kaup  and 
Newell  [14],  Karpman  and  Maslov  [13],  Knickerbocker  and  Newell  [19], 
Grimshaw  [3],  Kodama  and  Ablowitz  [17],  and  Smyth  [28]  show  that  a 
"shelf"  exists  there  and  thus  the  prediction  of  our  analysis  is 


consistent. 


We  now  return  to  the  discussion  of  the  leading  tail  of  a  right 


going  wave.  We  wish  to  solve  equation  (3.40)  subject  to  matching  the 
leading  order  singular  solution  (3.30).  Details  of  the  matching  will  be 
presented  in  Section  3.3  where  we  will  show  that  the  desired  matching 
solution  of  (3.40)  is  a  superposition  of  exponentials: 


*q(Z,t)  =/_ 


D(cc)exp[-ctZ 


3*(1V1)a2r]da, 


(3.41) 


where 


D(ci)  = 


(• 

|  0  exp[12a( 


x<'V1! 


V] 


a  <  0 
a  >  0  , 


and  D  is  determined  by  matching.  Equation  (3.41)  obviously  solves  the 
diffusion  equation.  The  substitution 


a  - 


-1/3. 

6X(lV'fy  6  ’ 


into  (3.41)  yields 


.  E-1/3bB  r 

=  6X(1)K(1!  ;o  exp[- 


6eAv1'K'1y 


/  2/3,, 

\-e  4? 


1 

+  2 


1/3h  .2 

e  bx5 


1  3X  . 


(3.42) 


a  convergent  integral  since  a  <  0.  It  is  an  exponential  integral  with  a 
cubic  phase,  an  incomplete  Airy  integral  of  the  type  which  would  arise 
in  the  analysis  of  the  formation  of  a  shock  wave  for  a  piston  (Lighthill 
[22],  although  no  details  are  presented). 

Using  Laplace's  method  for  the  asymptotics  of  exponential 
integrals  we  will  show  in  Section  3.5  that  (3.42)  is  the  desired 
matching  solution.  In  addition  we  will  show  that  the  critical  points  of 
the  exponential  phase  satisfy  the  fundamental  folding  incomplete 


quadratic  (3.18). 

We  point  out  here  that  if  we  solve  equation  (3.40)  as  an  initial 
value  problem  (which,  it  is  not),  with  the  initial  condition  given  by 
equation  (3.41)  with  x  =  0,  then  the  fundamental  solution  of  the 
diffusion  equation  yields  the  solution  we  have  obtained  by  matching. 

For  details  of  this  relationship  see  the  Appendix. 


3.3  Asymptotics  and  the  "Shock  Path” 

We  now  show  that  the  leading  order  solution  near  first  focusing 


=  *  exp[-K(1)(x-Av1  V15 


t)] 


matches  to  the  leading  order  singular  solution 
UQ  =  A  exp[-H]  . 


(3.43) 


(3.44) 


The  basic  matching  principle  is  to  express  each  solution  in  terms  of  the 
other  variables,  expand  for  small  e,  and  equate  the  resulting  expansions 
to  obtain  the  matching  conditions  (see  for  example  [15,  23]).  To  change 
variables  we  use  the  scaling  equations: 


x 


2/3 

£ 


(  t-t1  ) 


9 


and 

Z  =  e"2/3[X-X1-Cg(T-T1)-C(T-T1)Z] 

S  £1/3[x-x1-Cg(t-t1)-  C(t-t1 )2]  , 


from  Section  3.4.  Solving  for  T  and  X  we  have 
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Thus,  for  small  e  we  have 


H-i  [a'^  +  K(1  '(e2/^Z  *  e1/3C  T  .  e2/3Ct2) 
£  a 


XC1'K<i;3(e1/3x)l  , 


(3.45) 


and  from  (3.29) 


A  %  AU)|[b  +  2af(S)]e1/3tr1/2 


From  (3.42)  we  have 


(3.46) 


=  e"1//3d0/"exp[|-  hU)]d£  , 


(3.47) 


where 


TTTjrrr  > 

6Xv1/Kv1/ 


h(0  =  ~  a?3  +ib(T-T1)C2  -  [X-X1 -Cg(T-T1  )-C(T-T1  )2]?  . 

We  now  use  Laplace's  method  for  the  asymptotics  of  exponential  integrals 
on  (3.47).  The  critical  points  of  the  cubic  phase  are  determined  by 
solving 

h'Uc)  =  a£2  +  b(T-T1)^c  -  [X-XrCg(T-T1)-C(T-T1)2]  =  0  ,  (3.48) 

for  £,c.  Equation  (3.48)  is  the  fundamental  folding  incomplete 
quadratic  equation  (3.18)!  Using  (3.21)  and  solving  (3.48)  we  have 


5c  =  -2ir(T-V[1  *  ^  +  rf  s  ]’  • 

□ 


(3.49) 


Since  h(^)  is  a  cubic  with  a  negative  leading  coefficient,  its  local 
maximum  must  be  to  the  right  of  its  local  minimum  (when  they  exist  for 


V* _  ■*_  mf -  A A  if-a'V  n  n  *-*i «.  AWi  ^  *-*»-*■  “V  *-*i *- 


E 


1  -  /I  +  -r  S,  T  <  T  and  S  <  0 

=  -  ~(T-T,)  ,  bZ  1 

max  Zb  I  i  2 

1  -  /l  +  ^  S  ,  T  >  T  and  S  i 

'oL  1  43 


We  now  have 


,  l-1+3(1  +  S)  -  2(1  +  ■—  S)3/2 ,  T<T  and  S  £  0 

h(«  >  ■  ^2(t-t,)5  b  b 

max  2  1  I  ? 

I  “I  t  A  43  C  ^  oM  43  CN3/2  TNT  ^  c<-  ~b 

\-1+3(1  +  — ^S)  +  2(1  +  — 2s'  »  T>T  -jand  S±  — 


h'(£  )  =  0  , 

smax  ’ 


h"(r  )  =  b(T-T. ) 
^max  1 


✓l  +  -y  S  ,  T  <  T1  and  5  0 

- - b -  2 

-/  +  S  ,  T  >  T1  ,  and  5  -  ^. 

b 


In  Figure  3.8  we  have  extended  the  domain  of  h  to  illustrate  its  evolu¬ 
tion  throughout  the  neighborhood  of  first  focusing.  The  location  of 
5  and  the  magnitude  of  h(C  )  are  critical  in  the  determination  of 
the  asymptotics  of  using  Laplace's  method.  We  must  keep  in  mind 
that  £=0  is  a  cut-off  value  for  the  integral  in  4>g.  This  cut-off  gives 
rise  to  a  penumbral  caustic  in  contrast  to  a  more  common  cusped  caustic 
(for  further  discussion  of  the  classification  of  caustics  we  refer  the 
reader  ;o  [10]).  In  order  for  the  asymptotics  to  be  affected  by  £ 

m  3x 

we  must  have  ?max  >  0  and  b(^max)  —  h(0)  =  0,  otherwise  the  asymptotics 
will  be  completely  determined  by  the  endpoint  of  the  integral,  £  =  0. 
Laplace's  method,  for  fixed  5  >  0  and  T  <  T^  in  (3.47),  yields 


.  a/. 


Figure  3.8  Evolution  of  exponential  phase,  hU),  near  firs 
first  focusing  with  domain  of  h  extended  to 
negative  '  for  illustration  only.  Laplace's  method 
for  equation  (3.42)  uses  the  maximum  of  h(-;)  for  -  >  0. 


i 


40 


*n  '  s'1/3dC/"  expdthJC. .  )  .  V'U  )(C-C  )Z])d5 


9  1/3  H 

— - -  ]  1/2exo  NT 

h"(5  )J  P  Le  Cmax/J  * 

max 


(3.50) 


Now  using  (3.45),  (3.46)  and  (3.50)  in  equations  (3.43)  and  (3.44) 
yields  the  matching  condition 

D  =  A(0)[2e1/3dTT]"1/2  exp{l[e(  1 '-K(1 '(X  -X(  1  V1  ^  (3.51) 

e  I 

Having  thus  determined  5  we  can  now  obtain  the  asymptotics  of  ^  for 
T  >  T.  and  fixed  S  £-  b2/4a.  For  fixed  S  <  -3b_/(l6a),  5  *  £ 

I  l  iTlaX 

dominates  the  asymptotics  (see  Figure  3.8)  as  before  and  Laplace's 

2 

method  again  yields  (3.50).  However,  for  5  fixed  and  -3b  /(16a)  <  S  < 
2 

-b  /(4a),  the  endpoint  5=0,  dominates  the  asymptotics  (see  Figure 
3.8)  and  thus  using  Laplace's  method  for  an  endpoint  contribution  we 


find  from  (3.47)  that 


f>  'V 

u  Z  * 


(3.52) 


The  most  interesting  case  is  when  both  5=0  and  5  =  5max  >  0  contri¬ 
bute  to  the  asymptotics.  This  occurs  when  ^(5mgx)  =  h(0)  =  0  or 
2 

S  =  -3b  /(16a).  Laplace's  method  then  yields 


1/3 

d,  J  nr2e  d  1V2 

90  1  D[  h"(5  )] 

smax 


(3.53) 


This  clearly  demonstrates  that,  near  first  focusing,  there  is  a  transi¬ 
tion  region  where  the  asymptotic  solution  is  determined  both  by  the 

moving  core,  5=5  >0,  and  by  the  initial  conditions,  5=0.  This  in 

max 

fact  is  the  initiation  of  a  "shock  path".  The  solution  remains 
continuous,  but  there  is  an  abrupt  change  in  wave  number  or  a  wave 


number  shock.  Howard  and  Kopell  [10]  have  discussed  a  similar  concept 
for  slowly  varying  solutions  of  reaction-diffusion  equations.  The  shock 


path  near  first  focusing  (and  T>T^ )  is  determined  from 


Ss  =  [*s-XrCq(T-V-C(T-V2l/!r-Ii:2  =  -  Tfa  • 


the  points  where  the  asymptotic  are  determined  from  both  the  moving  core 
and  the  initial  conditions.  The  shock  path  near  first  focusing  is  thus 


Xs  ■  X1  *Cg(,-V  ♦  (C  -T?I;!T-V2  > 


(3.54) 


with  shock  speed 


2(C 


^(T-T  ) 
16a/U  V 


(3.55) 


As  we  entered  the  region  near  first  focusing  (T  <  T^ ) ,  we  matched 
the  leading  tail  solution  to  the  solution  near  first  focusing.  This 
matching  determined  the  solution  near  first  focusing.  We  must  match 
this  solution  to  the  solution  after  first  focusing  (T  >  T^ ) .  This 
matching  is  identical  to  the  matching  prior  to  first  focusing.  Thus  we 
see  that,  after  first  focusing,  the  moving  core  conLributes  a  single 
exponentially  decaying  tail  to  the  leading  tail  solution.  Since  in  this 
chapter  we  are  considering  the  case  with  no  crossing  of  characteristics 
due  to  the  initial  conditions,  it  is  easy  to  show  that  these  initial 
conditions  also  contribute  a  single  exponentially  decaying  tail  to  the 
leading  tail  solution.  Therefore,  for  each  space-time  poini  ahead  of 
the  core  where  a  single  characteristic  passes,  the  leading  order 
solution  is  a  single  exponentially  decaying  tail.  However,  inside  the 
envelope  of  crossing  characteristics  (see  Figure  3.3)  superposition 
yields  the  leading  order  solution: 
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u(x,t)  A  exp[-H  ]  +  Ar  exp[-Hr]  , 


(3.56) 


where  the  subscripts  "Jt"  and  "r"  denote  the  solutions  due  to  the  moving 
core  and  initial  conditions  respectively.  Equation  (3.56)  is  a  linear 
multi-phase  slowly  varying  tail  (see  [9,  24]).  The  asymptotic  solution 
for  the  leading  tail  is  a  single  tail  everywhere  except  at  or  near  where 


the  phases  in  (3.56)  are  equal,  or 


H  =  H 
l  r 

Applying  the  chain  rule  to  9  =  eH  we  have 


i!®.  -  k—  -  Q 
dT  ‘  aT 


(3.57) 


(3.59) 


Taking  X  =  X  (T)  to  be  the  points  where  the  two  phases  are  equal, 


we  have 


dXq  dXq 

K  _ _£  -  o  =  k  -  a  , 
l  dT  l  r  dT  r  ’ 


V  ar  [«] 


dT  '  K.  -  K 
*  r 


(3.59) 


This  is  the  well  known  Rankine-Hugoniot  relation  for  conservation  of 
waves,  equation  (2.24).  Using  (2.32),  equation  (3.59)  becomes 


VJ/'  c  7  o 

-3T2  =  +  K£Kr  +  K‘)  ,  (3.60) 

and  as  T  -*■  T.  we  see  that  (3.55)  and  (3.59)  do  match  since 
C  =  3Xv1/Kv  j  . 

g 

Figure  3.9  illustrates  the  asymptotic  solution  near  the  wave  number 


stock. 
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3.6  Constant  Coefficients 

We  now  illustrate  our  analysis  for  the  special  case  where  v,  X 
and  Y  are  constant  in  equation  (2.1).  Thus  (2.20)  becomes 


"dT”  =  "  3VW0  ’ 


wq(T)  =  Wg(0)exp[-  -jyT]  , 


(3.61) 


which  agrees  with  Karpman  and  Maslow  [13],  and  Ablowitz  and  Segur  [1]. 
The  leading  order  core  solution,  (2.17),  is  now  completely  determined, 
and  from  (2.13.b)  and  (3.61)  the  center  of  the  core  is  located  at 


vr:  *  4rt»o(o:  -  ”0(T)1  • 


(3.62) 


We  assume  X  >  0  and  thus  from  (2.17)  we  must  have  wQ(T)  >  0,  which  means 
the  core  moves  to  the  right. 

Now  turning  to  the  leading  tail  we  see  that  equation  (3.7) 

becomes 


dK  -  n 
ST  -  0  ’ 


(3.63) 


on  the  characteristics,  (3.8),  or 

K(X,T)  =  K[X  (c),?]  . 

0 

From  the  matching  condition,  (3.4),  to  leading  order  we  have 


wn  ^ '  i/2  2 

K(X,T)  =  [— ~ —  r/Z  exp[-^0  , 


XK4  =  wQ(5)  . 


(3.64) 


The  characteristics,  (3.0),  are  thus 


Figure  3.9  "Shock  path",  XS(T),  of  the  wave  number  shock  and 
asymptotic  solution  as  a  superposition  of  left  and  right 
tails  due  to  the  moving  core  and  initial  conditions. 


For  focusing  we  must  have 


3w'(r)(T  -r)  _  3wn(?)  +  X'  (0  =  0  , 


TfU)  =  ?  -  4y  (thus  Tf(C)  =  1  >  0) 


0.66) 


We  see  from  this  relation  that,  we  must  have  y  <  0  for  focusing  ahead  of 
the  core,  thus  the  core  must  accelerate!  First  focusing  occurs  at 


T  -  _  — 

1  "  Zy  ’ 


X1  =  " 


3wq(0/ 


From  (3.64)  we  have 


2AKKX  =  w'(c)Cx  , 
and  (3.65)  yields 


^C,  /  \  I  -ii  /- 


so  that 


3w6(0 


X  “  2  3w; 


f-5/-^w0u; 


(3.67) 


Using  (3.67)  in  (3.6)  yields 


A( X  ,  T )  =  A[X  (5), 5] 


2wq(c) 

3w^(5;(T-f)-2w0(5) 


exo[-y(T-c)] 


and  the  matching  condition,  (3.3),  then  yields 

12w  (?)  2w  '(g)  1/2 

A(X,T)  =  - - - T-rrrTTT  -s  exp[-Y(T-0]  .  (3.68) 

v  3w'  ( t,»  v l-t,  ,-2w  (  r, 


From  equations  (2.23)  we  have 


e(x,T)  =  k [ X— x  (?)]  _  n(r-s)  +  six  U),s]  , 

o  o 

on  the  characteristics,  or  from  (2.32)  and  matching  to  the  core 


6(X,T)  =  K[X-X0(C)  -  XlC(T-c)] 


(5.69) 


Thus,  from  (3.68)  and  (3.69),  the  leading  order  singular  solution  for 
the  leading  tail  is 


Ug  =  A  exp[-  . 


5.70) 


From  (3.48)  we  see  that  ?  satisfies  the  fundamental  foldinq 

max  y 

incomplete  quadratic  so  that  the  matching  solution  in  the  leading  tail 
(3.43),  becomes 


V  D[  ~"-2i"|-'b"(T-T~  )  ^/Z  exp[-K(1)(x-XK(1)Zt)]. 


(3.71) 


Matching  (3.70)  and  (3.71)  we  obtain  the  condition 
12w  (0)  ew  (0) 

B  =  — S - [J4-]172  , 

v  e'^TTd 


so  that. 


un  =  11  ]1/2  rT  exp[-^(-  e2/^Z  5  +  lE1/3btC2  +  -kc3)]dS 


V  £7T 


'0 


Cl'  T i ^ 

<p[-KU/(x-XKU/  t]. 


(3.72) 


Using  (3.23)  and  the  definitions  of  C^,  C,  a  and  b  we  find  from 


(3.65)  that 


Cg  =  3w0(0)  ’ 


C  =  0  , 

b  =  3wi ( 0 )  =  -4y  w„ ( 0 )  , 


Substituting  these  values  into  (3.54)  and  (3.55)  we  find  the  initial 
shock  path  and  speed  are 


CHAPTER  IV 


CUBIC  EOLDIMG/CUSPED  CAUSTICS 

4.1  Int roduct ion 

In  Chapter  III  we  considered  the  case  of  focusing  due  to  the 
moving  core  of  a  slowly  varying  solitary  wave  with  first  focusing 
occurring  along  the  characteristic  £  =  0  and  T^.(Q)  >  0.  We  found  for 
this  case  that  the  fundamental  folding  equation  was  quadratic  in  5  (see 
Figure  3.7),  and  gave  rise  to  a  penumbral  caustic  (see  Figure  3.4  and 
[20]).  We  now  turn  our  attention  to  the  other  two  possibilities:  Section 

4.2  gives  the  analysis  for  first  focusing  due  to  the  moving  core  along  a 
characteristic  >  0  and  T^(£,j)  =  0,and  Section  4.3  gives  the 
analysis  for  first  focusing  due  to  initial  conditions.  In  both  of  these 
cases  we  find  that  the  fundamental  folding  equation  is  cubic  in  the 
characteristic  variable  (see  Figure  4.2),  and  gives  rise  to  a  cusped 
caustic  (see  Figure  4.1). 

4.2  First  Focusing  along  ^  >  0  and  Tj.^)  =  0 

In  Section  3.2  we  found  that  the  characteristics  for  the  leading 
tail  due  to  the  moving  core  were  given  by 


X  =  F(T;£),  F(£;^)  =  X  (?),  £>0  , 


(4.1) 
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•  =  3AK  . 


For  focusing  we  have 


F?(Tf;C)  =  0  , 


Tf  =  Tf (? )  . 


(4.3) 


(4.4) 


(4.5) 


Thus  the  time,  ,  of  first  focusing  is  determined  by  minimizing  T^  for 
5^.0.  Chapter  III  analyzed  the  case  where  =  T^(0)  and  TL(0)  >  0. 

We  now  consider  the  case  where 


T1  =  W  * 


Tf(C.,)  =  0  . 


51  >  0  , 


Using  (4.7)  in  equation  (3.19)  we  find  that 


(4.6) 


(4.7) 


(4.8) 


so  that  the  fundamental  folding  equation,  (3.18),  must  include  cubic 


terms: 


6  FTTT(VV<T-V  =  FU^  ri  ;F1  '^T_T1  '(5-5,  ' 


.If  (t,;5,)(t-t,:2(Mi: 


+  2  FTSCT1  ’^1  '*■  T-ri  '  '  +  J  *^1 ' 


The  change  of  variable 


.V'.v.v.'.-;.-,'.-. j-.  .■ v .  ■  \ -\  -r.  \  ■*.  ■  <■. 


5 


€  "  Si  =  lT“Ti 


'f(S)  , 


(4.9) 


then  yields  the  similarity  variable 


S  =  [X-XrC  (T-T1)]/|T-T1  y!L  , 


and  the  fundamental  folding  cubic  is 


(4.10) 


af3(S)  +  bf(S)  =  S  . 


(4.11  ) 


For  this  case  we  have 


Cg  =  ft<t,!Ci:  , 


WWW  - 


*■  = 

By  differentiating  equation  (3.19)  with  respect  to  5,  we  find  that  a  <  0 
and  b  >  0.  From  equation  (4.11)  we  find  that  there  is  one  real  root,  and 
thus  from  (4.9)  one  characteristic  value  at  each  point  (X,T)  if 

Ta  '  vTa'  ?  U  ’ 

two  distinct  real  roots  (one  a  double  root),  and  thus  two  charac¬ 
teristic  values  at  each  point  (X,T)  if 

(2a/  +  '  -  0  ’ 

and  three  distinct  real  roots,  and  thus  three  characteristic  values  at 
each  point  (X,T)  if 

(2a'  +  v3a'  <  0  * 


Therefore,  using  (4.10),  the  cusped  caustic  is  given  by 


X  =  X1  +  C  ( T-T1 )  ±  2a[^|  (T-T.1 )  ]3/2  . 


(4.12) 


Figure  4.1  shows  this  cusped  caustic  and  the  number  of  characteristics 
in  each  region  due  to  the  moving  core,  while  Figure  4.2  gives  the 
evolution  of  the  number  of  characteristics  at  each  X  due  to  the  moving 


Proceeding  as  in  Section  3.3,  but  with  the  new  similarity 
variable,  (4.10),  we  find  that 

K  =  K(Tl5  &,)  +  i<T(T1 ;  c, )  (T-T1 )  +  it  (T.,  ,C1 )  f(S)  |  T-^  1 1/2+  . 


A  ^  S(?)|[b-3af2(S)](T-T1)|"1/2  . 

Using  these  expressions  in  (3.33)  we  then  find  that 


(4.13) 


B  |T-T1  |~5/2  f(S)  , 

and  continuing  to  higher  orders  shows  that  the  expansion  for  U  breaks 


down  when 


T-T1  =  0U1/2)  . 


Our  new  scales  are  thus 


x  =  e"1/2(T-T.)  =  e1/2(t-t.)  , 


Z  =  £"3/4tX-X1-Cg(T-T1)]  =  e1/4[x-xrC  (t-t.,)]. 
Now  seeking  a  solution  of  the  form  (3.38)  we  have 


(4.14) 


(4.15) 


r,(lUl)\  1/4, 


1/2  . 


t*v. 


Figure  4.1  Cubic  folding,  cusped  caustic,  and  the  number 
of  characteristics  at  each  X  due  to  the  moving  core. 


Thus  from  (2.21)  we  have 


The  substitution 


55 


into  (3.41)  now  yields 


i>0  =  e"1,/4d  D/“q  exp[—  h(^)]d  , 


(4.16) 


where 


— r-n-  f-n—  » 

6*v  1 'KU  ' 


h(c)  =  j  a(5-Cl)4  +  j  bd-T^Cc-^)2 
-  [X-XrCg(T-T1)](?-C1)  . 

The  integral,  (4.16),  with  quartic  phase,  is  similar  in  spirit  to  the 
oscillatory  Pearcey  integral  which  arises  for  cusped  caustics  in  diffrac¬ 
tion  (see  [9]).  Laplace's  method  for  the  asymptotics  of  this  exponential 
integral  shows  that  the  critical  points  of  the  quartic  phase  determined 


h'(Sc)  =  aUc-^r  +  b(T-T1)(Cc-C1)  -  [X-X^C  (T-^  )]  =  0  , 

satisfy  the  fundamental  folding  cubic  (4.11).  We  see  that  there  are 
one,  two,  or  three  real  distinct  critical  points  which  correspond  to  a 
maximum,  a  maximum  and  a  double  point  of  inflection,  or  two  relative 
maximums  and  a  relative  minimum  of  the  quartic  phase.  Figure  4.3 
illustrates  the  evolution  of  the  quartic  phase  in  the  region  near  first 
focusing  (the  cusp  of  the  caustic).  The  asymptotics  of  (4.16)  are 
determined  as  in  Chapter  III.  We  find  that  the  asymptotic  solution  is  a 
single,  slowly  varying,  exponentially  decaying  tail  everywhere  except  at 
or  near  the  points  where  h(£)  has  two  maximums  of  the  same  magnitude. 

The  points  where  this  occurs  yield  the  initiation  of  a  "shock  path"  as 


in  the  last  chapter,  where  the  asymptotic  solution  is  the  sum  of  two 
slowly  vary  tails  and  is  in  transition  from  one  to  the  other.  After  the 
wave  number  shock  forms,  the  shock  path  is  determined  by  the  inter¬ 
section  of  two  characteristics.  For  a  while  both  of  these  charac¬ 
teristics  initiate  from  the  moving  core  (5  >_Q).  However,  if  the 
characteristic  on  the  right  (5  <  5^)  tends  to  zero,  then  afterwards  the 
shock  path  is  determined  by  the  intersection  of  two  characteristics,  one 
from  the  moving  core  (5  >  5^)  and  one  from  the  initial  tail  (t  >  0). 

We  have  now  considered  two  of  three  possible  cases  for  focusing 
due  to  the  moving  core.  In  Chapter  III  we  showed  the  analysis  for  first 
focusing  along  the  characteristic  £  =  0,  with  TJ,(0)  >  0,  which  led  to 
a  fundamental  folding  incomplete  quadratic  and  a  penumbral  caustic.  We 
have  just  shown  the  analysis  for  first  focusing  along  the  characteristic 
£  r  C  >  0,  with  =  0,  which  led  to  a  fundamental  folding  cubic 

and  a  cusped  caustic.  The  final  possibility  for  focusing  due  to 
the  moving  core  is  where  first  focusing  occurs  along  the  characteristic 
5  =  0,  with  T£.(0)  =  0.  The  analysis  for  this  case  yields  a 
fundamental  folding  incomplete  cubic  (due  to  TJ.(0)  =  0,  see  Section 
3.3),  but  the  caustic  is  penumbral  because  we  have  5^.0  and  we  would 
need  C  <  0  to  generate  a  cusp.  Thus,  cutting  off  the  characeristic 
values  at  =  0  is  eauivalent,  in  this  case,  to  cutting  off  one  side  of 
the  cusp.  This  case  of  an  incomplete  cubic  can  be  easily  analyzed  by  our 
earlier  methods.  We  do  not  include  these  details  since  we  judge  this 
case  to  be  of  little  physical  interest. 


4.3  Focusing  Due  to  Initial  Conditions 


In  Section  3.2  we  found  that  the  characteristics  for  the  leading 


tail  due  to  initial  conditions  were  given  by 


X  =  G(T;;\  G(0; ; )  =  ',  c  >  0  , 
where,  as  before,  the  method  of  characteristics  yields 


ST  -  xK  ’ 


(4.18) 


dX  -  W 
ST  =  3  K  * 


From  equation  (4.18)  we  see  that  focusing  occurs  for 


G5(Tf;C)  =  0  , 


or  at  times 


T f  =  TfCO  . 


(4.19) 


(4.20) 


The  remaining  analysis  for  this  case  is  identical  to  that  in  Section  4.2 
for  the  previous  case  of  focusing.  The  only  differences  in  the  two 
cases  are  seen  by  comparing  Figures  4.1  and  4.2  with  Figures  4.4  and  4.5 
respectively,  which  show  that  the  characteristic  leading  to  first 
focusing,  £ ^ ,  eminates  from  the  initial  conditions  instead  of  the  moving 
core,  and  the  evolution  of  the  number  of  characteristics  at  each  X  is 
slightly  changed  since  all  of  the  characteristics  originate  at  time 


T  =  0. 


For  a  specialized  treatment  of  focusing  in  the  leading  tail  due 


to  initial  conditions  (where  slow  variations  leading  to  focusing  are  intro 
uced  strictly  through  slowly  varying  initial  conditions  and  not  by 
slowly  varying  coefficients)  the  reader  is  referred  to  [8]. 


APPENDIX 


1  1 

wflgl 


a 


Matching  Solution  of  the  Diffusion  Equation  Obtained  from  the 


Fundamental  Solution 


We  have  obtained  a  rather  unusual  solution,  (3.41),  of  the 


diffusion  equation,  (3.40).  Here,  we  will  show  that  there  are  initial 


conditions  which  will  yield  this  solution  from  the  fundamental  solution 


of  the  diffusion  equation: 


,(Z,  )  = 


«4tt(3X(1)K(1)t 


f  4>0(8,0) 


-(z-sr 

6XP  6  * 


(A.1) 


We  take  the  initial  condition,  4>q(  3,0),  from  the  matching  solution, 


equation  (3.41),  which  we  have  derived: 


4>Q(B,0)  =  O/gexpC-aS-iicr  ]da  , 


(A. 2) 


A(1)k(1)  2 

where  u  =  -12a( - r - ) L  . 


Substituting  (A. 2)  into  (A.1)  and  interchanging  the  order  of 


integration  we  obtain 


$g(Z,x)  =  5/gQ(Z,x)exp[-via  ]da, 


(A. 3) 


/ 4^(3 1  V(1)  )t 


-00  l  L  JJ 


-  ( Z- 3; 

exp  - TTTTTT —  ds 

4(3An  V1;)t 


(A. 4) 


However,  (A. 4)  is  the  fundamental  solution  of  the  following  initial 
value  problem: 

QT  =  3A(1)K(1)QZZ  , 


Q(Z  ,0)  =  exp[-aZ] . 


Thus,  we  solve  for  Q  and  obtain 


Q(Z,x)  =  exp[-az  +  3A(1  V1  )a2x]  . 


(A. 5) 


Substituting  (A. 5)  into  (A. 3)  then  yields  the  solution  of  the  diffusion 
equation  which  we  have  obtained  by  matching 


J>g  ( Z  ,  t)  =  D  J^expC-aZ  +  3Av^Kv^o^t  -ua^]da. 


(A. 6) 


BIBLIOGRAPHY 


[1] 

[2] 

[3]  R 

[4]  C 

[3]  R 

[6  J  R 

[7]  R 

[8]  R 

[9]  R 

[10]  L 

[11]  R 

[12]  T 

[13]  V 


.  0.  Ablowitz  and  H.  Segur  (1981),  Solitons  and  the  Inverse 
Scattering  Transform,  SIAM,  Philadelphia. 

.  0.  Ablowitz,  D.  0.  Kaup,  A.  C.  Newell  and  H.  Segur  (1974), 

The  inverse  scattering  transform-Fourier  analysis  for 
nonlinear  problems,  Stud.  Appl.  Math.,  53:  249-315. 

.  Courant  and  K.  0.  Friedrichs  (1948),  Supersonic  Flow  and 
Shock  Waves,  Interscience,  N.Y. 

.  S.  Gardner,  J.  M.  Greene,  M.  D.  Kruskal  and  R.  M.  Miura 

(1967),  Method  for  solving  the  Korteweg-deVries  equation, 
Phys.  Rev.  Lett.,  19:  1095-1097. 

.  Grimshaw  (1979),  Slowly  varying  solitary  waves.  I.  Korteweg- 
deVries  equation,  Proc.  Royal  Soc.  London  A,  368:  359-375. 

.  Grimshaw  (1978),  J.  Fluid  Mech. ,  86:  415-431. 

.  Grimshaw  (1970),  J.  Fluid  Mech. ,  42:  639-656. 

.  Haberman  and  D.  Allgaier  (1985),  Slowly  varying  solitary  wave 
tails:  focusing,  cusped  caustics,  wave  number  shocks  and 
birth  of  tails,  SIAM  J.  Appl.  Math.,  45:  919-927. 

.  Haberman  and  R.  Sun  (1985),  Nonlinear  cusped  caustics  for 
dispersive  waves,  Stud.  Appl.  Math.,  72:  1-37. 

.  N.  Howard  and  N.  Kopell  (1977),  Slowly  varying  waves  and  shock 
structures  in  reaction-diffusion  equations,  Stud.  Appl. 

Math. ,  56:  95-145. 

.  S.  Johnson  (1973),  On  an  asymptotic  solution  of  the  Korteweg- 
deVries  equation  with  slowly  varying  coefficients,  J.  Fluid 
Mech. ,  69:  813-824. 

.  Kakutani  (1971),  J.  Phys.  Soc.  Japan,  30:  272-276. 

.  I.  Karpman  and  F.  M.  Maslov  (1978),  Structure  of  tails 
produced  under  the  action  of  perturbations  on  solitons, 

Sov.  Phys.  JETP ,  48:  252-259. 


[14]  D 

[15]  3 

[16]  K 

[17]  K 

[18]  Y 

[19]  C 

[20]  Y 

[21]  S 

[22]  M 

[23]  A 

[24]  A 

[25]  K 

[26]  L 

[27]  3 

[28]  N 


.  J.  Kaup  and  A.  C.  Newell  '1978),  Solitons  as  particles, 
oscillators,  and  in  slowly  changing  media:  a  singular 
perturbation  theory,  Proc.  Roval  Soc.  London  A,  361: 

413-446. 

.  Kevorkian  and  3.  D.  Cole  (1981),  Perturbation  Methods  in 
Applied  Mathematics,  Springer-Verlag,  N.Y. 

.  Ko  and  H.  H.  Kuehl  (1980),  Energy  loss  of  a  Korteweg-deVr ies 
solitary  wave  in  a  slowly  varying  medium,  Phys.  Fluids, 

23:  834-836. 

.  Ko  and  H.  H.  Kuehl  (1978),  Korteweg-deVries  soliton  in  a 
slowly  varying  medium,  Phys.  Rev.  Lett.,  40:  233-236. 

.  Kodama  and  M.  3.  Ablowitz  (1981),  Perturbations  of  solitons 
and  solitary  waves,  Stud.  Appl.  Math.,  64:  225-245. 

.  3.  Knickerbocker  and  A.  C.  Newell  (1980),  Shelves  and  the 
Korteweg-deVries  equation,  3.  Fluid  Mech. ,  98:  803-818. 

.  A.  Kravtsov  and  Y.  I.  Orlov  (1983),  Caustics,  catastrophes, 
and  wave  fields,  Sov.  Phys.  Usp.  26:  1038-1058. 

.  Leibovich  and  G.  D.  Randall  (1973),  Amplification  and  decay  of 
long  nonlinear  waves,  3.  Fluid  Mech. ,  58:  481-493. 

.  3.  Lighthill  (1956),  Viscosity  effects  in  soundwaves  of  finite 
amplitude,  in  G.  K.  Batchelor  and  R.  M.  Davies,  eds. , 

Surveys  in  Mechanics,  250-351,  Cambridge  University  Press, 
Cambridge. 

.  H.  Nayfeh  (1981),  Introduction  to  Perturbation  Techniques, 
Wiley,  N.Y.  ~  “ 

.  C.  Newell  (1983),  The  history  of  the  soliton,  3.  Appl.  Mech., 
50:  1127-1138. 

.  Nishikawa  and  P.  K.  Kaw  (1975),  Phys.  Lett.  A,  50:  455-456. 

.  A.  Ostrovsky  and  E.  N.  Pelinovsky  (1970),  Izv.  Acad.  Sci. 

USSR  atmos.  oceanic  Phys.,  6:  552-555  (English  translation). 

.  S.  Russell  (1844),  Report  on  waves,  Report  of  the  14th 

Meeting  of  the  British  Association  for  the  Advancement  of 
Science ,  3ohn  Wiley,  London.  ~ 

.  F.  Smyth  (1985),  Formation  of  solitons  for  a  perturbated 
Korteweg-deVries  equation,  preprint. 


[29]  G.  B.  Whitham  (1974),  Linear  and  Nonlinear  Waves,  Wiley,  N.Y. 


